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ABSTRACT 

We present the serendipitous discovery of 195 extragalactic globular clusters (GCs) in one of the 
deepest optical images ever obtained, a 126 orbit Hubble Space Telescope (HST) Advanced Camera 
for Surveys (ACS) imaging study of the nearby Galactic GC NGC 6397. The distant GCs are all 
found surrounding a bright elliptical galaxy in the field, and are among the faintest objects detected 
in the image, with magnitudes 26 < F814W < 30. We measure the redshift of the parent elliptical 
galaxy, using the Gemini Multi-Object Spectrograph (GMOS) on Gemini South, to be z = 0.089 
(375 Mpc). This galaxy, and its associated clusters, therefore ranks as one of the most distant such 
systems discovered to date. The measured light from these clusters was emitted 1.2 Gyr ago (the 
lookback time) and therefore the optical properties hold clues for understanding the evolution of GCs 
over the past Gyr. We measure the color function of the bright GCs and find that both a blue 
and red population exist, and that the colors of each sub-population are redder than GCs in local 
elliptical galaxies of comparable luminosity. For the blue clusters, the observed color difference from 
z = 0.089 to today is only slightly larger than predictions from stellar evolution (e.g., changes in the 
luminosity and color of the main-sequence turnoff and the morphology of the horizontal branch) . A 
larger color difference is found in the red clusters, possibly suggesting that they are very metal-rich 
and/or significantly younger than 12 Gyr. 

Subject headings: galaxies: elliptical and lenticular, cD - galaxies: star clusters - globular clusters: 
general - stars: evolution 



1. INTRODUCTION 

The formation of globular star clusters (GCs) is closely 
linked to the evolutionary histories of their parent galax- 
ies. Major epochs of star formation in isolated galaxies 
as well as interactions with other gala xies lead to sub- 
sequent spikes in the n umber of GC s (ISchweizerl 119871: 
Ashman fc Zepl fl998t lHarrisI 12001 iBrodie & Straderl 
2006f ). Indeed, it appears that the processes involved in 
galaxy formation imprint a signature on the clusters as 
reflected by the similar properties of the host galaxy and 
its associated clusters (e.g., age and metallicity) . The 
prevalence of a large number of GCs around all mas- 
sive galaxies, as well as their optically bright luminosi- 
ties (—10 < My < —5), makes them excellent tracers of 
galaxy formation processes. 

To date, GCs have been detected in hundreds of nearby 
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galaxies (see Table 1 in Brodie & Strader 2006 for those 
systems with accurate photometry). However, direct 
age measurements from resolved main-sequence turnoff 
fitting are only possible for the nearest clusters in the 
Milky Way, LMC, SMC, and M31 (e.g., see the recent 
study of SKHB 312 by Brown et al. 2004). Properties 
of more distant systems are typically constrained with 
integrated light photometry and, where possible, spec- 
troscopy. Such work has demonstrated that the luminos- 
ity and color distributions of GCs in most galaxies are 
remarkably similar (e.g., Harris & Racine 1979), with 
most populations exhibiting a bimodal color distribution 
reflecting old metal-poor and metal- rich subpopulations. 
Whether these populations reflect GCs that formed at 
two different times (e.g., Forbes et al. 1997) is still de- 
bated. Unfortunately, directly probing the properties of 
these clusters at earlier times is a difficult task. For ex- 
ample, at the distance of the nearby Coma cluster, a 
typical GC with My — — 7.5 has an apparent magnitude 
of V = 27.5. 

In HST Cycle 13, our team (PI: Richer - GO- 
10424) was granted 126 orbits of ACS time in a single 
pointing to explore the white dwarf cooling sequence, 
faint main-sequence, and space m otion of the Milky 
Way core-collapsed GC NGC 6397 (iRicher et alJ 120061 : 
lHansenet all 120071 : iKalirai et ai1l2007al) . The resulting 
data set represents one of the deepest optical images ever 
obtained in astronomy, with a 50% (25%) completeness 
limit of F8UW = 28 (F8UW = 28.75) and a photo- 
metric error of a = 0.3 at F8UW = 28.0 and .F606IU 
= 29.25. In this Letter, we present the discovery of a 
distant, background galaxy in this data set that har- 
bors its own GC population. We measure the redshift 
of the galaxy to be z — 0.089, and therefore this system 
resides ~375 Mpc from the Milky Way at a lookback 
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Fig. 1 — A 3.4' X 1.7' section of the 126 orbit HST exposure of NGC 6397 (this is approximately one-half of the total ACS field of 
view). The strong gradient across the image from top- left to bottom-right reflects the decreasing radial density of NGC 6397 stars (the 
cluster center is located ~6.5 arcminutes to the N-W of our field). The inset shows a closer look (8" X 8" — 15.9 X 15.9 kpc at the distance 
of the elliptical galaxy) of the large elliptical galaxy at image coordinates (x, y) = (789.9, 750.9) pixels. The point source enhancement 
surrounding the elliptical galaxy is very high and suggests that these objects are GCs themselves, orbiting this galaxy. Over the area of 
the inset, only one cluster white dwarf is expected to contaminate our sample. Five of the GCs are highlighted with small circles in the 
inset. Higher resolution version of this figure is available at http://www.ucolick.org/>~jkalirai/DistantGlobs/ Image credit: NASA, ESA, 
H. Richer (UBC), J. Kalirai (UCSC). 

time of 1.2 Gyr. The color function of the GCs in this 
galaxy is found to show differences when compared to 
local samples suggesting the possible detection of evolu- 
tionary changes in GC properties over the past Gyr. 

2. IMAGING AND SPECTROSCOPIC OBSERVATIONS 

The HST/ ACS imaging observations of NGC 6397 
from GO-10424 have been described in several papers 
(e.g., Richer et al. 2006; Hansen et al. 2007; Kalirai et al. 
2007a; Anderson et al. 2008). To summarize, we ob- 
tained 252 images in F814W and 126 in F606W-", for 
a total integration time of 126 orbits in a single point- 
ing centered at aj2ooo = 17:41:03, 5j2ooo = —53:44:21. 
Sources detected at the same position on 90 of the 252 
.F814W images were shown to be a 3<r detection (from 
simulations) and cleaned to remove false detections in 
the wings of bright stars, those caused by intersecting 
diffraction spikes, and objects morphologically inconsis- 
tent with stars (e.g., galaxies and cosmic rays). PSF 
p hotometry on these sou rces was performed as described 
in lAnderson et all (J2008D . 

Our expectation was that the final cleaned catalog of 
^8000 objects from this analysis represented stars on the 
image. However, a spatial plot of the faint-blue sources 
indicated a clustering of 244 objects centered on a bright 
F814W = 16.7 elliptical galaxy at pixel location x = 
789.9, y = 750.9 (see Figure □}. We reduced the ini- 



tial catalogue by 6 objects that look extended, another 
5 objects that are very bright and red, and 38 objects 
that were only detected in one filter and therefore for 
which we have no color information. This final sample of 
195 objects must be GCs themselves, in orbit around the 
background galaxy. The clusters are among the faintest 
objects in our catalogue, ranging in brightness from 26 
< F814W < 30. Given the ACS pixel scale and distance 
to the galaxy (see below), these clusters are all located 
within a radial distance of 1.7 - 20 kpc (0.85" - 10") 
from the center of this galaxy. The incompleteness in 
our sample increases rapidly as we approach the center 
of the galaxy and therefore we preferentially detect the 
bluer clusters. We can roughly estimate the masses of the 
clusters assuming they are old (12 Gyr) and metal-poor 
([F e/H] = —1 . 5). F or a Kroupa initial-mass function, 
the iMarastonl (|2005f ) models indicate that the GCs with 
F814W = 27. 5 are ~2 x 10 6 M© whereas the brightest 
globulars in this galaxy are ~10 7 M@. 

To measure the distance to the galaxy accurately, we 
first obtained a ground based image of the field in excel- 
lent seeing conditions and low air-mass with the Gemini 
Multi-Object Spectrograph (GMOS) on Gemini South 
(Program ID: GS-2006A-DD-16, executed on 18 Aug. 
2006). Despite the confusion from NGC 6397 stars, 
the elliptical galaxy was easily delineated at aj2ooo — 
17:41:12.2, 5j2ooo = —53:43:16.5 in a 120 s exposure with 
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the g' filter. Next, we obtained a low resolution spectrum 
of the galaxy with a 1" longslit and the B600 grating (R 
= 1600) centered at 5200 A. The position angle was set 
to 325 degrees to avoid contamination in the slit from 
nearby stars. Four 900 s target exposures, as well as ob- 
servations of a flux standard, a Cu/Ar arc, and a flat field 
were obtained. The data were reduced using the Gemini 
IRAF Package, Version 1.8. 

The final reduced spectrum of the elliptical galaxy 
in the observed wavelength frame of reference is shown 
in Figure [2] (top). The arrows mark several well- 
defined spectral features along with their rest-frame 
wavelengths.. A cross-correlation with a template spec- 
trum of an old, metal-rich single stellar population yields 
a best-fit redshift of z = 0.0894. Assuming standard 
concordance cosmology with Q m — 0.3, fl\ = 0.7, and 
h = 0.7, this redshift corresponds to a proper radial dis- 
tance of 375 Mpc and a lookback time of 1.2 Gyr. The 
luminosity distance of 409 Mpc gives a distance modulus 
of m — M — 38.06. While this is not the most distant 
galaxy in which GCs have been detected (that currently 
belongs to Abell 1689 at z = 0.183; Mieske et al. 2004), 
the photometry in this galaxy is considerably deeper than 
in any previous study of a galaxy at a comparable dis- 
tance. 

3. RESULTS - PHOTOMETRIC PROPERTIES OF 
DISTANT GLOBULAR CLUSTERS 

To probe the properties of these distant GCs, we com- 
pare their color distribution to that of globulars in lo- 
cal elliptical galaxies. These galaxies nearly always have 
a bimodal distribution of GC colors, with peaks near 
V — I ~ 0.95 and 1.18. There is a mild dependence of 
the peak GC colors with host galaxy luminosity (Strader, 
Brodie, & Forbes 2004; Peng et al. 2006), and so we cor- 
rect these values to V — I ~ 0.93 and 1.16 for a galaxy 
with M v = -20.5 (assuming V - I = 1.2). We note 
that our measurement of the integrated brightness of the 
galaxy is only roughly correct given the large number of 
foreground contaminating stars in the outer parts. This 
does not effect our results as a 0.5 magnitude difference 
in the galaxy luminosity would translate to a <0.01 mean 
magnitude offset in the expected peak color of the GCs. 

In Figure [2] (middle) we present the color magnitude 
diagram of 196 GCs. This includes the foreground glob- 
ular NGC 6397 where we resolve all of the individual 
stars (small points) as well as the 195 extragalactic ob- 
jects (larger points below the white dwarf cooling sequence 
of NGC 6391). A large fraction of these distant GCs 
(~60%) have very faint magnitudes and lie below the 
50% completeness limit (dashed line, F814W^ = 28). A 
total of 46 objects, or ~64 after accounting for incom- 
pleteness, have F81AW < 27.5. The foreground redden- 
ing in t he direction of NGC 6397 is E(F606W - F8UW) 
= 0.18 ([Hansen et al.ll2007n . giving A F814W = 0.33 for a 
standard reddening curve. Using the distance modulus 
in § 2 and this reddening, the expected peak of the ap- 
proximately lognormal GC luminosity function (GCLF) 
is at F814W ~ 30.0. It follows that we have only ob- 
served the tip of the GCLF in this galaxy. Assuming a 
standard a = 1.3 for the GCLF (e.g., Kundu & Whit- 
more 2001), down to F8UW = 27.5 only -2.7% of the 
total GC population is visible (assuming a Gaussian form 
for the entire GCLF). Even though this estimate is very 
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Fig. 2. — Top - A Gemini GMOS spectrum of the elliptical galaxy 
shown in Figure [T] reveals its redshift to be z = 0.0894, indicating 
a luminosity distance of 409 Mpc. Several absorption features and 
their rest-frame wavelengths are illustrated. Middle - The proper 
motion selected color-magnitude diagram of the foreground cluster 
NGC 6397 shows a beautiful main-sequence extending from the 
turnoff down to the lowest mass stars that burn hydrogen (red 
points). The rich white dwarf cooling sequence of the cluster is also 
visible. The larger dots in the faint-blue part of the CMD represent 
the 195 extragalactic GCs found in the vicinity of the large elliptical 
galaxy in Figure [l] Bottom — Histogram and density estimate for 
the colors of 46 GCs with F814W < 27.5. The majority of GCs 
belong to the blue subpopulation, with a tail towards red objects. 



rough, it appears that there must be more than a thou- 
sand GCs in this galaxy. This puts the galaxy squarely 
in the regime of massive ellipticals, and we can safely use 
the peak GC colors of local ellipticals as a reference. 

In Figure [2] (bottom) we show a color histogram of GC 
candidates with F814W < 27.5 and a density estimate 
overplotted. There is a broad peak of blue objects with 
F606W - F8UW between -0.8 and 1.2, and a tail to 
redder values. Mixture modeling of the color distribu- 
tion with Nmix (see discussion in Strader et al. 2006) 
finds that the preferred fit is two Gaussians centered at 
F606W-F8UW = 1.03 ± 0.02 and 1.29 ± 0.05. The er- 
rors are derived through bootstrapping. While increasing 
photometric errors tend to broaden the color distribution 
as one goes to fainter magnitudes, performing the same 
fitting to a limit of F81AW = 28 yields essentially the 
same result: a preference for two populations of GCs, 
with peaks at F606W - F8UW = 1.03 and 1.31. The 
small number of red GCs makes that peak value uncer- 
tain, but the derived blue peak appears consistent with 
a visual estimate from Figure [2] (bottom). 

We now compare these observed peak values to those 
expected from a scaling of local GC systems in a galaxy 
with same luminosity as our host elliptical. To convert 
from V — I to F606W — F814 W, we use the empirical 
relation of lHarris et al.l (|2007f ) , derived from photome- 
try of NGC 2419. This relation is F606W - F8UW = 
0.707(y - I) + 0.033(V - I) 2 . To these values we add 
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Galactic reddening of E(F606W - F8UW) = 0.18±0.03. 
Finally, we add a /c-correction derived from Bruzual & 
Chariot (2003) stellar population models. This depends 
slightly on the (unknown) metallicity of the stellar pop- 
ulation, but assuming values of [Fe/H] ~ —1.2 and —0.2 
for the two subpopulations gives fc-corrections of 0.088 
and 0.063, respectively. If our assumed metallicities 
are incorrect, the effect on the F606W — F81AW col- 
ors is minor: a 0.1 dex change in [Fe/H] corresponds to 
about 0.01 mag, and 0.2 dex is a reasonable assump- 
tion for the uncertainty in the metallicity. Combining 
these corrections together gives predicted peak values 
of F606W - F8UW = 0.95 ± 0.03 and 1.11 ± 0.04 
disregarding evolution in the stellar population. We 
also computed the predicted peak colors assuming the 
iSirianni et al.l (2005) transformations between Vega mag- 
nitudes and ACS filters and found the same results, 
within the small error bars. 

The observed change in the peak colors of blue and red 
GCs from z = 0.089 to today is 0.08 ± 0.04 and 0.18 ± 
0.0 6 magnitudes. W hat is the expected evolution? Both 
the lMarastor] (|2005l) and lBruzual fc Charlotl (|2003l ) mod- 
els predict that the red metal-rich GCs should be about 
0.01 mag redder in F606W — F814W at the present day, 
and that there should be essentially no (< 0.01 mag) 
evolution in the blue metal-poor GCs. Therefore, to first 
order, our results suggest that more evolution has oc- 
curred than is predicted by the models, especially for 
the red clusters 8 . We note that some GC systems in the 
local universe show the peculiar property of a correla- 
tion between blue GC magnitude a nd color, such that 
the more luminous GCs are redder ( Harris et al.l 120061 : 
IStrader et aDl2006l : iMieske et alJl2006f ). Since we sample 
only the brightest GCs in this distant galaxy, the mean 
color of the blue GCs derived above may be redder than 
the typical object, therefore bringing the predicted and 
observed values closer to the expected difference. We 
therefore conclude that the behavior of the blue GCs 
is essentially consistent with simple stellar population 
models, while the red GCs are substantially redder at a 
lookback time of ~ 1.2 Gyr than expected. 

4. DISCUSSION 

The most obvious candidate for any change in blue 
GC colors is the horizontal branch (HB). At old ages, its 
morphology is highly sensitive to a number of variables, 
including age. At younger ages, GCs should typically 
have red HBs, which transition to blue HBs as the stellar 
population ages. For old systems, the quantitative evolu- 
tion in a relatively red color such as F606W — F81AW is 
expected to be small. Our results confirm this, suggest- 
ing that blue GCs have become bluer b y just 0.08 ± 0.04 
magnitudes over the past 1.2 Gyr. The iMarastonl (|20Q5l ) 
models for a 12-13 Gyr metal-poor ([Z/H] = —1.35) sin- 
gle stellar population predict a difference of only 0.02 
mag between an object with a blue and a red HB, con- 
sistent with our result at 1.5o\ Observations of Galactic 
GCs show no correlation between residuals oiV—I colors 
from a fiducial V — I to [Fe/H] relation and two different 
parameterizations of HB morphology, suggesting that — 
at least for the range of HBs observed in the Galaxy — the 

8 Based on the tightness of NGC 6397's main-sequence, the dif- 
ferential reddening along this line of sight is negligible. 



HB has little effect on the V — I color (Smith & Strader 
2007). 

For the red GCs, our mixture modeling fit of the color 
distribution indicates a peak at F606W - F81AW = 1.29 
± 0.05, 0.18 ± 0.06 magnitudes redder than the local 
sample. Although our data set contains few red GCs, 
the observed peak is found to be at the same color ir- 
respective of the magnitude limit adopted (i.e., pushing 
the sample to F8UW < 28). Over the past 1.2 Gyr, the 
expected evolution from models of a single stellar pop- 
ulation is much smaller than this, as noted above. One 
possible way to reconcile some of these observed differ- 
ences is to consider whether a fraction of these clusters 
are significantly younger than 12 Gyr, and more metal- 
rich. In this case, the predicted color difference result- 
ing from standard stellar evolution between a redshift of 
Z = 0.089 and now would be larger than that calculated 
above. 

An additional effect may explain the observed color 
difference if the meta llicities of the cluster s are in fact 
supersolar. Recently, iKalirai et al.l (|2007bl ) have shown 
that the remnant population of white dwarfs in the 
old (8 Gyr), supersolar metallicity Galactic star clus- 
ter NGC 6791 are undermassive when compared to more 
metal-poor systems. These stars likely formed through 
a unique channel involving enhanced mass loss of the 
progenitor star on the red giant branch. Such evolution 
naturally explains both the absence of stars near the tip 
of the red giant branch in this cluster and the population 
of the extreme horizontal branch of the system (see also 
Castellani & Castellani 1993). Although the morphology 
of the horizontal branch does not heavily affect the color 
evolution of a population in our filters, the depletion of 
stars near the tip of the red giant branch can be more 
important and lead to an increasingly bluer population 
as the system ages. Of course, this mechanism would 
be even more efficient if the clusters are older since the 
envelope mass of the evolving stars is lower. However, 
this me chanism does require the clusters to have [Fe/H] 
> +0.2 (jKilic. Stanek. fc Pinsonneaultll2007h . and we see 
few systems nearby that are this metal-rich. 

Finally, we note that these objects are far enough away 
that we cannot rule out a contaminating population of 
ultra-compact dwarfs (UCDs). The properties of these 
potentially transition objects (between GCs and dwarf 
galaxies) vary widely among galaxies. For example, in 
NGC 3311, they are are seen as very red objects more 
luminous than the old and red GC population (Wehner 
& Harris 2007). However, in NGC 1407, the UCDs have 
larger sizes than globulars but have colors that fall be- 
tween the blue and red subpopulations (Harris et al. 
2006). We note that UCDs are typically found in massive 
galaxy clusters or "groups", and it is difficult to assess 
the environment near our elliptical galaxy given the pres- 
ence of NGC 6397 stars out to several arcminutes (1' = 
120 kpc). 

5. CONCLUSIONS 

We have presented a study of the photometric prop- 
erties of one of the most distant sample of GCs yet dis- 
covered, 195 objects at d = 400 Mpc. These clusters are 
all found around a bright elliptical galaxy that is itself 
located behind one of the nearest Galactic GCs to the 
Sun, NGC 6397. Given their large distance, the light 
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from these clusters has taken f.2 Gyr to reach us, and 
therefore these objects offer us a rare chance to probe the 
evolution of GCs in the recent past. By comparing the 
color function of these clusters to nearby GCs, we find 
differences that suggest the mean colors of GCs have be- 
come bluer in the past Gyr. This evolutionary change 
can be explained assuming standard stellar evolution for 
the blue clusters. However, we find differences that sug- 
gest the mean colors of the red GCs were significantly 
redder 1.2 Gyr ago as compared to samples in local giant 
ellipticals, much redder than evolutionary models pre- 
dict. 
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